Symmetrized models of last passage percolation and non-intersecting lattice paths 



(N 



Peter J. Forrester and Eric M. Rains 

*Department of Mathematics and Statistics, University of Melbourne, 
Victoria 3010, Australia ; 

^ Department of Mathematics, University of California, Davis, CA 95616, USA 

K" 

' It has been shown that the last passage time in certain symmetrized models of directed per- 

, eolation can be written in terms of averages over random matrices from the classical groups 

U{1), Sp{2l) and 0{l). We present a theory of such results based on non-intersecting lattice 
paths, and integration techniques familiar from the theory of random matrices. Detailed 
derivations of probabilities relating to two further symmetrizations are also given. 

^ '. 1 Introduction 

■ There are a number of striking results linking models of stochastic processes to to random matrix theory 

(for a recent work of this type see [1]; for reviews see [71 [H]). As an easy to explain example, consider 



the following last passage percolation problem due to Hammersley. In the unit square mark in points 
uniformly at random according to a Poisson rate with intensity A (thus the probability the square contains 
N points is equal to -^e^ ). Join points by straight line segments with the requirement that the segments 
have positive slope and form a continuous path, and extend this path to begin at (0, 0) and finish at (1, 1). 
1^ , Define the length of the extended path as the number of points it contains, and denote by — 1^{X) the 

lO ! stochastic variable specifying the maximum of the lengths of all possible extended paths (see Figure [T]). 

' Then it is known from the work of Gessel [TT] and Rains [12] (see also [S]) that 

cn 

Vn ■ Pr(/° < = /TTe^^^'^A (1.1) 

O ■ ^ ^ =1 

where the average is with respect to the eigenvalue probability density function (p.d.f.) of random 
matrices chosen uniformly at random from the group U{1). The latter has the explicit form 



o 



Jwi n Ie"^--e"^r (1.2) 



<i<fe<; 

The formula (|l.ip is central to the proof by Baik, Deift and Johansson [T] giving the limiting scaled 
distribution of the longest increasing subsequence of a random permutation. 

In a substantial work Baik and Rains have considered symmetrized generalizations of the Hammer- 
sley process, and have shown that the cumulative probability for the analogue of the stochastic variable 
P can be written as an average over the classical groups Sp{2l), 0(1) (for two particular symmetries), 
or U{21), U{1) U{1) (for two other symmetries). Moreover, random matrix formulas were also given [3] 
for lattice generalizations of these processes (in the case of the original Hammersley process the lattice 
generalization is referred to as the Johansson model [13]). The proofs of f2] make sophisticated use of 
symmetric function theory, and have independent interest in that discipline. On the other hand it is 
reasonable to suggest that many researchers interested in directed percolation and growth processes will 
lack the necessary background in symmetric function theory to fully appreciate these proofs. This has 
motivated us to seek alternative derivations which to a large extent avoid the heavy use of symmetric 
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Figure 1: Eight points in the unit square, and the extended directed paths of maximum length. 
Since the number of segments in these paths equals four, here /° = 3. 

function theory called upon in [2] . Instead our proofs make use of a non- intersecting path picture of the 
Johansson growth model [TH H] , an extension of this picture to a last passage percolation model with 
Bernoulli random variables, together with techniques familiar from the theory of random matrices (the 
applicability of such techniques have been foreshadowed in Section 6 of [5]). Another essential ingredient 
from [2\ is the use of bijections from the theory of Young tableaux (see e.g. [IS])- Here our presentation 
differs from that of [2 only in that we give more detail. 

Our task then is to derive, from a non-intersecting paths picture, formulas known from [2] giving 
cumulative probabilities of a suitable last passage percolation variable for generalizations of the Ham- 
mersely process in terms of random matrix averages. As already mentioned, one such generalization is 
the Johansson model. In Section 2 its definition is recalled, as is its formulation in the non-intersecting 
path picture. Formulas from the theory of non-intersecting paths are then used to derive the analogue 
of (jl.ip . An analogous discussion of a variant of this model, involving Bernoulli rather than geometric 
random variables, is given in Section 3. The four symmetrized versions of the Johansson model are then 
treated separately in each of the subsequent four sections. 

2 Johansson model and polynuclear growth 

We begin by defining the last passage percolation model introduced by Johansson [T3]. Consider the 
upper right quadrant square lattice : i,j £ Z+}. Associate with each lattice site a random 

non- negative integer variable Xi.j chosen from the geometric distribution with parameter athj so that 



Denote by (1, l)u/rh(n,n) a sequence of lattice paths starting at (1, 1) and finishing at {n,n) with each 
lattice point successively connected by edges which are either directed upwards or horizontally to the 
right. One defines the last passage time say as the maximum of the sum of the integer variables 
associated with these lattice points. Thus 



fc) = (1 - a^hj){a,hjf. 



(2.1) 




— max 




(2.2) 



(l,l)u/rh(n,n) 
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Figure 2: RSK correspondence in the polynuclear growth model picture between a weighted 
non-negative integer matrix and a pair of weighted non-intersecting lattice paths. 

According to [2], for given parameters {ui}, {bj} the cumulative distribution can be written as a random 
matrix average according to 

n n I 

Pr(£° < = n (1 - «^^^-)( n 11(1 + a,e-'"^'')(l + b.e''^'')) (2.3) 
jj=i j=i k=i ^ ' 

where the average over U{1) refers to the probability density function p.2|l . We seek a derivation of 
(|2.3p . and analogous formulas from 2i for symmetrized versions of the Johansson model, within a non- 
intersecting paths representation of the model. 

Let us first revise how non-intersecting paths relate to the Johansson model [HI [9] . This is done via 
a geometrical construction, equivalent to the Robinson-Schensted-Knuth (RSK) correspondence from the 
theory of Young tableaux |10j . which gives a bijection between non- negative integer matrices and pairs of 
non-intersecting lattice paths. Furthermore, the maximum displacement of the top-most of these paths 
is equal to L^, with the profile of this path also specifying the height profile in a statistical mechanical 
model referred to as the polynuclear growth model. With regard to the latter, the entries Xij of an n x n 
non-negative integer matrix X = [xi,j]ij^i^,,,^„ (for convenience rows are labelled from the bottom) now 
represent the heights of columns of unit width centred about x — j — i which occur at time t = i + j — 1 
(in labelling the matrix in terms of x and t it is convenient to first rotate it 45° anti-clockwise). The 
columns arc to be placed on top of the profile formed by earlier nucleation events and their growth. The 
right boundary of the column corresponding to Xij is to be weighted a^'"^ while the left boundary is 
to be weighted b^''\ and these weights are to be multiplied together with any existing weights along 
the same vertical segment of the profile. During each time interval the existing profile or profiles are 
required to grow one unit to the left and one unit to the right, with any resulting overlap, together with 
the corresponding portion of the weights, recorded on a profile with base one unit below. In this way 
a bijection between n x n integer matrices with each entry Xi_j weighted aibj, and a pair of weighted 
non- intersecting paths is obtained. A particular example is given in Figure O 

The bijection generates at most n non-intersecting paths. These paths start one unit apart in the 
y-direction (at y = 0, (ri — I)). In the x-direction one member of the pair of paths starts at 
X = — (2n — 3/2) and may go up in arbitrary integer amounts of a unit step at x = — (2n -f 1/2 — 2j) 
for j — 1, ... ,71 with each step weighted by bj, or to the right in steps of two units (these steps are 
unweighted). The other member starts at a; = (2n— 3/2) and may go up in arbitrary integer amounts of 
a unit step at x — (2n + 1/2 — 2i) for i = 1, . . . , rt with each step weighted Oi, or to the left in steps of 
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two units, the latter being unweighted. Furthermore the second member is constrained so that it joins 
with the first member at a; = 0. Both members are equivalent to what will be termed u/rh (up/ right 
horizontal) non-intersecting lattice paths. By definition such paths are defined on the square lattice and 
start at a; = 0, one unit apart in the y-direction (at y = 0, . . . (tt. — 1)), and finish at x = n — 1, 
with y-coordinates iii — {I — 1) {I = 1, . . . ,n) where /^i > /i2 > • • • > > 0. The path starting at 
y = — (Z — 1) is referred to as the level-/ path. Each path may move either up or to the right along the 
edges of the lattice, with each up step at x = j — 1 weighted qj. Define the weight of a configuration of 
u/rh lattice paths as the product of all the step weights. Then it is well known (see e.g. [22]) that with 
:= (/ii, /i2, . . . , 11 n) (because of the orderings of the /i;, /i forms a partition) 

(weight of the paths) = s^(gi, . . . , g„) (2.4) 

u / rh. paths 
displacements fj. 

where is the Schur polynomial. 

We remark that u/rh non- intersecting lattice paths are equivalent to semi-standard tableaux (num- 
bered diagram of a partition A such that the numbers weakly increase along rows and strictly decrease 
down columns). Thus with A; (j) denoting the number of vertical steps in the level-Z path at x = j — 1, the 
Ith row of the tableaux is of length X]J=i ^lU) ^^i^ is numbered by Xi{j) lots of j's (j = 1, . . . , n in 

order). An explicit example is given in Figure [31 Consequently if there are n lattice paths the numbering 
is from the set {1, . . . ,n} which is referred to as the content of the tableau. With A = (Ai, . . . , A„) 
denoting the partition formed from the length of the rows, the tableaux is said to have shape A. 

For future reference we note that with /x;(n,j) denoting the displacement of the level-/ path at 
X — ~{2n + 1/2 — 2j) as resulting from the RSK correspondence and fii{i,n) equal to the displacement 
of the level-/ path at a; = 2n + 1/2 — 2i, we have 

n n 

J2iMnJ)-f^iin,j-l)) = ^x.j (2.5) 

1=1 i=l 
n n 

^{fii{i,n)~ fii{i-l,n)) = ^x^J. (2.6) 

1=1 j=i 

Also for future reference we make note of quantities generalizing which are related to the maximum 
displacements of the level-/ paths for each / = 1, 2, . . . . For this let (rd*)' denote the set of / disjoint 
(here meaning connecting no common lattice sites) rd* lattice paths, the latter defined as either a single 
point, or points connected by segments formed out of arbitrary positive integer multiples of steps to the 
right and steps up in the square lattice 1 < j < n. Generalizing the definition (j2.2p by defining 

L°«=max^x..„ (2.7) 

(rd*)' 



a theorem of Greene [12] gives 



and thus in particular jl5j 



I 



L°^'^ = E Z^"- (2-8) 



m— 1 



,.r = L°W:=L°. (2.9) 



It follows from the above discussion that with the entries of the n x n matrix X weighted according 
to (|2.ip , the probability that X maps under the RSK correspondence to a pair of non- intersecting paths 
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Figure 3: Drawn in heavy lines on the square lattice is a family of two u/rh lattice paths with 
up segments allowed at x = 1,2,3,4, while drawn in heavy dashed lines are the corresponding 
four dual lattice paths. Also given is the semi-standard tableau encoding the two u/rh lattice 
paths. 

with maximum displacements /i is given by |14| 

n 

J]^ (1 - aibj)Sf^{ai, . . . ,a„)sp(6i, . . . ,6„). (2-10) 
The equality (12. 9p between fii and then implies the formula [Tl] 

n 

Pr(L°<0 - Y[{l-a,bj) s^{ai,...,an)s^{bi,...,bn). (2.11) 

= l fi:fj,i<l 

Thus we must now show that the sum over Schur functions in (|2.11|) can be written as the average over 
U{1) in (|2.3p . Moreover we want to achieve this task within the framework of non- intersecting paths. 

An important notion for this purpose is that of the dual non-intersecting lattice paths associated with 
a set of u/rh paths. The dual lattice paths connect points on the lattice {(n — 1/2, m), n € Z>o, m G Z}. 
Points are connected by segments which are directed either right horizontal (rh) or diagonally down (dd) , 
with a dd segment bisecting every u segment of the u/rh lattice path. The dd segments are connected 
by rh segments to form dual lattice paths starting at x ~ —1/2 and finishing at x = n — 1/2 in the x- 
direction, while in the y-direction these paths start at y — k for fc = 1, 2, . . . , /zi where /zi is the maximum 
displacement of the level-1 path (see Figure [3]). In terms of tableaux, the dual lattice paths correspond 
to reading down columns instead of across rows. The mapping carries over to weighted paths by simply 
weighting the dd segment in the dual path by the value of the u segment it bisects in the original u/rh 
path. 

The key feature for our purpose is that the constraint on the u/rh lattice paths having maximum 
displacement less than or equal to / translates in the dual path picture as constraining the number 
of paths to be less than or equal to L As only dd segments are weighted we can take the number 
of paths as being exactly equal to I. Furthermore we can regard a pair of rh/dd lattice paths, each 
with the same end points and each containing I paths but weighted from {bj} and {a^} respectively, as 
a single set of non-intersecting paths. In the latter the second member of the pair is reflected about 
X — n — 1/2 so that its final positions are at (2n — 1/2, k) for k — 1,2, ... ,Z at it consists of right 
horizontal and up diagonal segments, the latter being weighted by aj according to them passing through 
X = 2rt + 1/2 — j. With the initial and final y coordinates generalized to y^\ . . . , J/;^'^'' and yi,. . . ,yi 
respectively, let G2n{yf \ ■ ■ ■ .. .,yi) denote the total weight of all such paths. It then follows 



from the dual lattice paths interpretation of the sum in (|2.1ip that 



Pr(L° <l)=l[{l- a^fij)G2n{{yf = -{j - {Vj = -(.? - (2.12) 

Because each path in the family is directed, the weight G2n of all paths in the family is given in terms of 
the weight of a single path in the family, g2n{y''^^', y) say. Thus the well known Linstrom-Gessel-Viennot 
theorem (see e.g. [TB]) gives 



G2n {yf^ , ■ • ■ , yf ; yi , . . . , y/ ) = det \ g2n (y)"' ; yk) 
Furthermore it is easy to see that 



,(0).. 



,(0). 



j,k=i,...,i 



(2.13) 



1 I-2TT N 

92n{y^°^;y) = ^ I + a,e-^<>){l + b,e^'>)e~^(y~y'°'^Ue. (2.14) 

Substituting (|2.14p in (|2.13[) , then substituting the result in (|2.12p and recalling the general formula [23] 



det 



27r 



j,A:— l,...,n 



U{n) 



(2.15) 



we see that (12.31) is reclaimed. 



3 Polynuclear growth with BernoulU random variables 

A variant of the last passage percolation model revised in the previous section is to replace (|2.ip by the 
Bernoulli distribution 

Pr(x,, = fc) = -^^i^. A: = 0,1. (3.1) 
Let X = [xij] i=i,...,rrv be an array of such variables. One specifies the corresponding last passage time by 

j=l,...,Tl 

J^m.n—max Xi,j, (3.2) 

(i' ,j')ebottom to top 
u/rd paths 

where the sum is over all u/rd paths in the rectangle 1 < i' < m, I < j' < n from the bottom row (row 
1) to the top row (row m). The segments of the path join entries successively to the north or north-east 
in the array. 

Underlying this model is the dual RSK correspondence 10 . To our knowledge this has not previously 
been related to a layered growth model formed out of non-intcrsccting paths. Here such a relationship 
will be presented. 

Again, the entries Xij of the array X are regarded as recording nucleation events. However, unlike 
the situation with the RSK correspondence itself, the entries of X are not first rotated 45° before being 
associated with positions and times. Rather x^.j = 1 denotes a nucleation event (a unit square) which is 
positioned above the segment x — j — ltox — j, and on top of earlier nucleation events and their growth 
(this is in common with the polynuclear growth model of the previous section) . These nucleation events 
occur at successive times < = 1, 2, . . . , n, with the positions recorded by I's in the corresponding rows of 
X. Thus to begin, at i = 1 the nucleation events are read off from the first row of X and marked on the 
line y — 0. As t t + 1, the existing profile(s) is to grow one unit to the right (but not the left) until 
it joins up with the neighbouring nucleation event on the right. If there is no such right neighbour, and 
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y=-2 



Figure 4: Mapping from a 0,1 matrix to a pair of non-intersecting lattice paths. 

this nucleation event is yet to grow (i.e. recorded in the previous time step), it is to grow to a: = n + 1 
and have its shape modified by removing the upper triangular half of its final square. If it has right edge 
at a; > n + 1, it is to grow one unit to the right. Also, the meeting of all nucleation events in going from 
t i-^ t + I are to be recorded on the line y = —t as new nucleation events with left edge at the positions 
of the meetings. 

This procedure is to stop after time m + I along y = —{I — 1) {I = 1, . . . ,to + 1), this being the 
maximum time for which new nucleation events can be created and then grow once. The layers of profiles 
which are so formed are of the form u/rh (up/ right horizontal) non-intersecting paths from a; = to 
X = n — I, and Id/lh (left diagonal/ left horizontal) non-intersecting paths from x = n + in to x = n + 1 
(see Figure HI). 

According to (|2.4p and the surrounding text the total weight of all non-intersecting u/rh paths initially 
equally spaced at y — 0, . . . , — (n — 1) along x — 0, finishing at y = /ij — (j — 1), (j/ = 1, . . . , n) along 
X = n — 1, with up steps at x = j — 1 weighted bj is given by s^(6i, . . . , 5„). Further, the dual paths of 
Figure[3]are precisely the non-intersecting Id/lh paths initially equally spaced at y = 0, . . . , — (n— 1) along 
X = n + m, finishing at ?/ = /ij — (j — 1) {j — 1, . . . , n) along a; = n -I- 1, which make up the second family 
in the growth process. With the possible up steps (each of unit length) atx^n + m + l — i weighted by 
Qi, the total weight of the paths is s^i{ai, . . . , &„), where /i' denotes the partition corresponding to the 
transpose of the diagram of fi. 

It follows from this that with an n x m array of O's and I's formed according to (|3.ip . the probability 
that it maps, under the dual RSK correspondence, to the above specified nested growth profile with 
maximum displacement /j, is given by 



m n 



YiYii^ + aibj) ^s^/(ai,.... 



am)s^(6i, . . . , 



bn) 



(3.3) 



(cf. (|2.10p ). Analogous to (|2.1ip we have the normalization condition 



l[l[{l + a,b,)-'Y.s^,{a,,..., 



am)Sf,{bi, . . . , bn) = 1. 



(3.4) 



i=ij=i fi 
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It follows immediately from (|3.3p that 

m n 

Pr(C,„ </) = [][](l + a,&j)-i J2 Sf,'{ai,...,am)s^{hi,...,bn). (3.5) 

From Figure [3] and related text, we know that the Schur polynomial s^i can be interpreted in terms 
of non-intersecting u/rh paths by reading along rows, while can be interpreted in terms of u/rd paths 
by reading down columns. In both cases the number of paths is equal to /zi, which according to (|3.5p is 
no bigger than I. Further, the left set of paths consists of m steps, and the right set consists of n steps. 
Let G*^„^{P'^;P'^), := (Z - 1, Z - 2, . . . , 0) denote the total weight of such paths. Analogous to ((21^ 
we have that 

m n 

Pr(i?„i„ < I) ^X{'Y{{l + a.h,)-^Gl^^(P^-}S^'>). (3.6) 

But according to the Linstrom-Gessel-Viennot theorem 

Gl^^P^.P^) = det[g:,„(/f ;4°))],,,.i,..., (3.7) 

where gi* ^{x,y) is the weight of a single path of the prescribed type starting at x and finishing at y. 
The latter can readily be seen to be given by 

9Ui^;y)^j- / l[{l + a,e^<^)l[{l-b,e~^<^r'e-^'^y'^U9 (3.8) 

j = l k=l 

Substituting (|3.8p in (|3.7p . making use of (|2.15p . and substituting in p.6p we reclaim the expression for 
Pr(i^„ < /) as a random matrix average [2], 

m n ml n I _ ^ 

p^Liin < I) - n 11(1 +«»^^)"'((n 11(1 (n 11(1 - 6,e-^^))" ) . (3.9) 

i=lj=l j=lk=l j=lk=l ' 

4 Matrices symmetric about the anti-diagonal 

With our convention of numbering rows from the bottom, the term anti-diagonal used here is what is 
conventionally referred to as the diagonal of the matrix. Under the RSK correspondence in the non- 
intersecting paths formulation, matrices symmetric about the anti-diagonal give a bijection with pairs of 
non- intersecting u/rd lattice paths in which only one member of the pair is independent. 

Regarding this point, reflect the entries of a general n x n non- negative integer X about the anti- 
diagonal to form the matrix = [xn+i-j.n+i-i]i,j=i,...,7i- We see from the definition (|2.7p that the 
quantities in^'^ are invariant under this transformation, and thus according to (|2.8p so then too are the 
path displacements /i;. Furthermore, it follows from (|2.5p that 



'^(f^finj) - fJ.fin,j - 1)) = ^Xr,+l-j,^ = ^{ni{n+l~ j,n) ~ ^i{n~j,n)) 

1=1 i=l 1=1 

n n n 

- (i - !,«)) = ^Xj^ri+i-i = ^{ni{n,n + I - j) - ni{n,n- j)). (4.1) 



1=1 1=1 1=1 



These equations respectively tell us that the total number of up steps at x = — 2n— i + 2j [x = 2n + ^ — 2i) 
in the paths corresponding to X^ is equal to the number of up steps at x = 2j — | ^ — 2i + |) 
in the paths corresponding to X. Furthermore, ii X = X^ an algorithm can be presented (due to 
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Schiitzenburger; see e.g. [2^) which aUows paths from a; = 2n — 1/2 to x = 1/2 to be constructed out 
of the paths from x = —2n + 1/2 to a; — —1/2. This permits a bijection between weighted matrices 
symmetric about the diagonal and a single set of u/rh lattice paths provided the weighting of Xi,j in 
the former is equal to (aia„+i_j)^''J (and thus 6n+i-i = Oj)- For example, the matrix of Figure [5] is 
symmetric about the anti-diagonal. With n = 2, setting bn+i-i = ai {i — 1, . . . ,n) we see that the weight 
of steps at a; = — 2n — 1/2 + 2i is equal to the weight of up steps at a; = 2i — 1/2 (i = 1, . . . , n). 

With the constraint X — X^, to obtain a probabilistic setting we should weight only the sites 
i < n + I — j . To achieve this without affecting the weights of the pairs of paths we simply square the 
weights at the sites i < n + 1 — j, and set the weights at sites i>n+l— jto unity. With a.; — ^Jqi^ we 
therefore choose 

Pr(a;ij = A:) = (1 - qiqn+\-j){<ii<ln\\-jt , i < n + I - 

Pr(a;,,„+i_, = fc) = (1 - q,)q^ (4.2) 

which corresponds to weighting the vertical segments of the single u/rh paths by gi, 52, . . . , <?n from right 
to left. According to (|2.4p with the maximum displacement of the level-/ path denoted by fii, the total 
weight of such paths is given by s^((;i, . . . , g„). We conclude that under the RSK mapping, with the 
independent entries chosen according to ()4.2|) . the probability a non-negative integer matrix symmetric 
about the anti-diagonal maps to a pair of u/rh paths with final displacement fi is equal to 

n 

n^^"**) n (1 ^ ft'jjOsMi'?!' • • ■ (4-3) 

The specification f|4.2p can be generalized, allowing for a generalization of (|4.3p . For this one recalls 
[15] that the RSK correspondence has the property that for [xij] symmetric about the anti-diagonal 

n 

#{xi,„+i_i : Xi^n+i^i odd} = : fij odd} = ^{~''^y~^ fJ-'j (4-4) 

where /i^- denotes the displacement of the level j conjugate path or equivalently the length of the jth 
column in the corresponding tableau (recall Figure [3]). Hence if we generalize the second probability in 
(021) to read 

Pr{x,,n+l-^ = k)= ^-^—^P" (4.5) 
1 + pg^ 



then we have that (14. 3p generalizes to 

II n (l-q.q,)f^"^^'-''''''^s,iq,,...,q^. 



1=1 ' ■ l<j<j<n 



Writing 



:= max Xi 



(l,l)u/rh(„,™) 



it follows from this that [2] 



Pr(L,^<0 = nrTj^ n (l-*9,) E (4-6) 

i=l ' l<i<j<n ^:Ml<i 

Using symmetric function theory, Baik and Rains [2J show that the sum in (|4.6p can be written as a 
random matrix average involving the classical group Sp{2l). Matrices from this subgroup of U{21) have 
their eigenvalues in complex conjugate pairs {e^*^^ }j=i^....;, with < 9j < tt {j — 1, . . . , I). Here we will 
give a derivation in keeping with integration techniques from random matrix theory. 
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Proposition 1. Consider the eigenvalues with angles < 9j < tt {j — . . . ,1) of matrices from Sp{2l). 
Let { )sp{2i) denote an average with respect to the corresponding eigenvalue p-d.f, 



1 1 



(4.7) 



i<j<fe<; 



One has 



f'(^-^^"=nTT|- n (i-«,)(n(^r37E^nii 

i^l ^ l<i<j<n k^l ' ^ ' j^l 

n I n 



i0k\2 



Sp{2l) 



l<i<j<n 



k=lj=l 



Sp(2/) 



(4.8) 



(4.9) 



Proof The maximum possible height in the growth model relating to (|3.5p is m, implying the so called 
dual Cauchy identity 



J]^ + a^bj) = ^ s^'(ai, . . . ,a„i)s^(6i, . . . ,6„). 

i—1 j — 1 fii<.7n 



(4.10) 



Renaming the parameters, it follows from this that 



fe=i j=i 



fi:fii<.2l 



(4.11) 



Substituting (|4.1ip in (|4.8p . substituting the result in (|4.6p . equating coefficients of Sp,{qi, . . . ,qn) and 
writing ^' ~ p shows that (|4.8p is equivalent to the matrix integral formula 



|l_^e-*e,|2^pie ,e ,...,e ,e j/^^^^n " ^ 



fe=i 



Sp(2l) 



(4.12) 



Now, the meaning of the matrix integral is an integral over 0^ e [0, tt], weighted by (j4.7p . Noting that the 
integrand is unchanged by 0j i— s- —9i, and making use of the determinant formula for Schur polynomials 



s\{xi, . . .,xn) 



det[g; 



JV-fc+Afci 



\j,k=l,...,N 



det[qf-%,k=i,...,N 



(4.13) 



shows that the matrix integral is equal to 

I 



1 



A;=l 



/3e- 



-jefc|2 



det 



g-«Sj(p2i-fc + l+fc-l) 



(4.14) 



3 = 1, 

k = l,...,2l 



The structure of the integral (|4.14p is very common in random matrix theory [18[ [6]. It can be written 
as a Pfaffian, giving that (|4.14p is equal to 



2' 



iS{pj~j-Pk+k) _ p~iO(pj-j-pk+k) 



.2-K |l-/3e-^^|2 
which after evaluating the integral therein reduces to 

/3"'Pf [sgn(fc - j)/3l«-'"'+'^-J'l 



)de 



j,k=l:---,2l 



j,k=l,...,2l 



(4.15) 



(4.16) 



10 



This Pfaffian is special case Xj = pj ~ j, f{xj) = in the general formula [8] 



Pf 



\f(x,)J 



sgn(a;j - Xk) 



n 



.f{xQi2j-l)) 



(4.17) 



where the permutation Q is such that 

XQ{2j-i) > Xq(2j), Q{2j) > Q{2j - 1) (j = 

and thus evaluating to the r.h.s. of (|4.12p . 

In regards to (|4.9p . use of an appropriate modification of (|4.1ip shows that this is equivalent to the 
matrix integral formula 



Sp{2l) 



To derive this, use of (|4.7p and (|4.13p shows that the matrix integral is equal to 
1 



(4.18) 



rPf 



fe=i 



I3e- 



gidj (P21-fc + 2 + fc-l) 
g-idj (P2!-fc + 2 + fc-l) 
'P2!-fc + 2 + fe— 1 



k = l,...,2l+l 



Ai2l+l)xi2l+l) + 



j=l,...,2l + l 



where 



^(2/+l)x(2/+l) 



1 



fe=i,...,2;+i 



27r |l-/3e- 



iS{Pj-j-pk+k) _ -i9{pj-j-pk+k} 



)de 



jM=l,...,2l + l 



Here the second equality follows from standard integration methods of random matrix theory, in the same 
way that (|4.15p follows from (|4.14p . Computing the integral reduces this to 



^-(/+i)pf 



[sgn(A: - ])p\P^-P'+''-%,k=i,...,2i+i [(i''=+^'+^-%=i, 



,21 



^-(iP'+^'+^-%^l,...,2l 

But this is precisely the same as (I4.16P with I ^ l + l, P21+2 — 0, and so reduces to the r.h.s. of (|4.18p . □ 



5 Matrices symmetric about the diagonal 

According to the rules of the polynuclear growth model, if a non-negative integer matrix X — 
maps to a pair of u/rh non-intersecting lattice paths (^1,^2) of the same final displacement, then the 
transposed matrix — maps to the pair of u/rh non-intersecting paths (P2,^'i)- Hence 

the Robinson-Schensted-Knuth correspondence when applied to symmetric matrices X — X"^ gives a 
bijection with a single set of u/rh lattice paths, since then we must have Pi ^ P2- To obtain a bijection 
between weighted symmetric matrices and a weighted set of u/rh lattice paths, and furthermore to obtain 
a probabilistic setting, we weight only the entries i < j, with the value of x^.j for i > j fixed by symmetry. 
Arguing then as in the derivation of (|4.3p we see with 

Pr(a;,j = fc) = (1 - q^qJ){q^qjf , i < j Pr(a;,,, = fc) = (1 - qi)q^ (5.1) 
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the probability the symmetric matrix X maps to a set of u/rh paths with final displacements ^ is equal 

to [13] 

n 

11(1-9') n {^-q^<ll)sMl.■■■,<ln)■ (5.2) 

As with (|4.3p this can be generalized to the case that the diagonal entries are chosen according to 

Pr(x,,, = fc)-(l-ag,0(ag^)^ (5.3) 
Thus recalling |15[ [H] that in the Robinson-Schensted-Knuth correspondence for symmetric matrices 

n n 

E^" = E(-i)'"V,' (5.4) 

with the generalization (|5.3p , (|5.2p should correspondingly be generalized to read [2] 

n 

\{{l-aq,) W {l~m,)a^'^-^(-^'>'"^^Sf,{qi,...,qn). (5.5) 

2=1 l<i<j<n 



Writing 



Lf^ := max ^ 



(l,l)u/rh(„,„) 



and noting that 

n I 

E(^l)"'^^Mj = #(columns of odd length in yu) = ^ /"'fcinod 2 

J = l /c=l 

where / = /ii, it follows from (|5.5p that 

n 

Pr(L^ < = - n (l-*9.) E a^^-^^^'"°'''sM(9i>---:9n). (5.6) 

Starting with this formula, Baik and Rains [5] proved the following analogue of Proposition [1] involving 
now a random matrix average involving the classical group 0{l). Matrices from this subgroup of U{1) 
form two disjoint components, 0~^{l) and 0~{l), distinguished by the determinant equalling +1 or — 1 
respectively. The complex eigenvalues occur in complex conjugate pairs, and there is a real eigenvalue 
z — —1 for matrices in 0^(1) with I odd, a real eigenvalue eigenvalue z = 1 for matrices in 0^{l) with I 
even, and two real eigenvalues z — ±1 for matrices in {I) with / even. 

Proposition 2. Consider the eigenvalues with angles < 9j < tt, {j = I, . . . , I) of matrices from 0{l). 
Define 

{ ■ )o{i) = 2 ' ^"^^^'^ ^ ^ ' 

where ( • )o+(i) denotes an average with respect to the eigenvalue p.d.f. for random matrices from the 
classical group 0^{l), 

,V22^-ia/2V n |e--e-N^|l-e^(«^-«^)p, ^ even (5.7) 
i<i<*;<'/2 
. ('-i)/2 

^(;-i)/22/-i(a_i)/2)l 11 1^ I 11 1^ e I |i e | , i oaa, 

" '' ] = 1 l<J<fe<(i-l)/2 

(5.8) 
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and (•)o-(/) denotes an average with respect to the eigenvalue p.d.f. for random matrices from the classical 
group (l), 

(5.9) 
(5.10) 

We have 

n n 

Pr(^^ < - 11(1 - "9') n {i-mj)(det{li + aU)ll{h+q,U)) (5.11) 

Proof. Use of the dual Cauchy identity (liTTU)) in ([g?TT|) and comparison with (j5.6p shows that (j5.1ip 
is equivalent to the matrix integral evaluation 

(det(li+aU)sJU)) ^a^5=iP.mod2 /5-^2) 

where Sp(C/) denotes the Schur polynomial as a function of all the eigenvalues of U. 

Consider first the I even case, I i-^ 21, and consider separately the components 0^{2l) of 0{2l). 
Substituting the eigenvalue p.d.f. for 0+(2Z) (|5.7p . and proceeding as in the derivation of (|4.16p . which 
involves use of (|4.13p and reduction to a Pfafhan, shows 

(det(l2i + aU)sp{U))ueo+{2i) = 2^''Pi[ajkhk=i,...,2i (5.13) 

where 

O'jk = ((1 + a^)'5(p^-i)-(pfc-fc),odd + 2a(5(p^_j)_(p^_fc),cvon)sgn(fc - j) 

= (i(l + a)2-i(l-a)2(-l)(«-^)-('"=-'=))sgn(fc-j). (5.14) 

The task is therefore to compute the Pfaffian of the matrix with these entries. For this one uses the 
identity [H] 

Pf(A + S)= (-l)^^«^'"'^'/'Pfs(^)Pfs(S) (5.15) 

SC{1,2,...,2Z} 
IS I even 

where Pf5(yl) denotes the Pfaffian of A restricted to rows and columns specified by the index set 5*, and 
similarly Pf5(B). With 



A 



.2 

and noting 



i(l + a)2sgn(fc-j)l . , S= f-i(l-a)2(-l)(«-^")-('"'-^)sgn(fc-j) 



j,k=l,...,2l 



2 



j,k=l,...,2l 



Pf[sgn(fc - j)] = 1, Pf[aj- fe(-l)(«-^')-('"=-'=)] = (-l)^(''^-^')pf[aj- fe] 

we see that 

PfsA = 2-1^1/2(1 + a f\, Pf^B = (-2)-l^l/2(l - a)2'-|s|(_i)E,«ft-j. 
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It thus follows from (|5.15p that 

1 , x-^ /I + a\ l-si /I - a\2i-|s| „ 

2i-'Pf(A + i?) = 2 (4^) (^) (-1)^^-^"^ 



(5.16) 



SC{1,2,...,2Z} 



Now, in general 



J2 x\S\y2l-\S\^_^^^^^sP. 



SC{l,2,...,2i} 
I SI even 

_ 1 

~ 2 



E 

SC{1,2,...,2;} 
21 



(_y)2;-|5|(_i)E,eS«) 



SC{1,2,...,2;} 

2/ 

Using this result to evaluate (|5.16p and substituting in (|5.13|) gives the matrix integral evaluation 



(5.17) 



We turn now to the corresponding formula for the average over O {21). The analogue of (|5.13p in 
this case is ^ 

(det(l2/ + aU)sp{U))u<,o-{2i) = ^^^r^[C]Pf K.fe + C&j-,fckfc=i,....2/ (5.18) 

where aj^k is as in (|5.13p while bj^k = {—1)P''^'^ — (— l)''^^^ , and [C] denotes the coefficient of (. Observing 
that 

[b,k] - uw^ - wu^, u = [l]j=i^...,2;, w = [(-l)''^-^],=i,...,2/ (5.19) 

shows that \bjk\ has rank 2. It follows that the Pfaffian in ()5.18p is linear in C, and so the r.h.s. of ()5.18p 
can be rewritten 

With 7, Ci, ^2 arbitrary non-zero constants, the structure of \hjk\ and use of elementary row and column 
operations verifies that this in turn can be rewritten 



(5.20) 



(l-g^) 1 
2'-i C1C2 



Pf 



[a-jk] Ciw C2U 
-(iw^ 7 
-C2W^ -7 



7Pf [a 



(5.21) 



Setting (i ~ i(l — a)^, C2 = (1 + c^)^, adding one half of the final row/column to the second 
last row/column, and subtracting the second last row/column from the final row column, then setting 
7 = (1 + a^) allows (|5.2ip to be rewritten as 



Pf[ajfc 



2(/+l)x2(/+l) 



P21 + 1—P2L+2—O 



(l + a^)Pf[ajfc]2/x2i 



l-a2 

Comparing (|5.13p and (|5.17p tells us that 

Substituting in (j5.22p and simplifying implies the matrix integral evaluation 

{dct{Ui+aU)sp{U))ueo-(2i) = a^" _ ^E?Li(P.+i)mod2 



(5.22) 



(5.23) 
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Adding this to (|5.17p verifies (|5.12p in the case / even. 

Similar working suffices in the I odd case, Z i-^ 2/ + 1. In regards to the average over 0^{2l + 1), 
making use of the exphcit form of the p.d.f. (|5.9p . the determinant form of the Schur polynomial (|4.13p . 
and integration techniques from random matrix theory, one obtains the Pfafhan formula 



(det(l2/+i +aC/)sp(C/))c/eo+(2;+i) - -^^^Pf 



[0'jk]j,k=l,...,2l + l [l]j=l,...,2/+l 
— [i-]k=l,...,2l + l 



(5.24) 



where is specified by (|5.14p . This Pfaffian can in fact be evaluated by making use of the Pfaffian 
evaluation implied by the equality of (|5.13p and (|5.17p . To see this, multiply the final row and column 
of the matrix in (|5.24p by (1 + a)^, and balance the equation by dividing by a prrefactor of (1 + a)^ on 
the r.h.s.. Next subtract the 2nd last row from the final row, and 2nd last column from the final column. 
Finally, write in the 2nd last entry of the final row and column (1 + a)^ = (1 + a^) + 2a. This shows 
that (|5.24p is equal to 

^77T^Pf[A2;+2)x(2Z+2) +B'] (5.25) 

where ^(2/+2)x(2/+2) ■= [0'jk]j.k=i....,2i+2 while B' has all entries zero except for the second last entry of 
the final column, which is 2a, and the second last entry of the final row, which is —2a. Making use of 
(j5.15p shows that this in turn is equal to 

fPf ^(2i+2)x(2i+2) 



2'(1 + a) 



+ 2aPfA2,x2i)- (5-26) 

P2;+2— P2i + 1 



But the value of Pf A21X21 for general I is known from the equality between (|5.13p and (|5.17p . so we find 
that (|07|) reduces to 

Q,E?LV P,mod2 ^ Q,E ■LY(P.+l)mod2 ^5 27) 

thus giving the evaluation of the random matrix average in (|5.24p . 

For the average over 0~{2l + 1), we note that a change of variables 0j 1-^ tt — thetaj {j = 1, . . . ,1) 
shows 

{detil2i+i + aU).Sp{U))ueo-(2i+i) = (-l)!"! (det(l2i+i - aC/)sp(t/))c/eo+(2i+i) • (5-28) 
Substituting (|5.27p with a t-^ —a for the average on the r.h.s. shows that this is equal to 

a^f=i Pjm°d2 _ ■LY(p.+i)m°d2_ ^5 29) 

Finally, taking the arithmetic mean of (|5.27p and (|5.29p . we obtain the sought evaluation (|5.12p with 
I 1-^21 + 1. □ 



6 Matrices symmetric about both the diagonal and anti-diagonal 

Let the 2n x 2n matrix X = [xij]ij=i^,,,^2n have non-negative integer entries, and label the rows from the 
bottom. Suppose furthermore the entries are symmetric with respect to reflections in both the diagonal 
{xij = Xj^i, i > j) and anti-diagonal (x^j = Xi^2n+i-j, i > 2n + 1 — j) . Because X is symmetric about 
the diagonal, the RSK correspondence maps X to a pair of identical non-intersecting u/rh lattice paths 
(P, P) say. On the other hand X being symmetric about the anti-diagonal implies X maps to the lattice 
path pair {P^ , P), where P^ is the Schiitzenberger dual of P. Consequently in this case X maps under 
the RSK correspondence to a single set of at most 2n u/rh non-intersecting lattice paths P with the 
special property that P = P^ . 
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We will suppose furthermore that the entries on the anti-diagonal are constrained to be even. Then 
according to (j4.4p all final displacements of the paths must be even. A partition with parts 2Xi so 
each part is even will be written 2A. 

The independent elements of X can be chosen to be Xi^j with z < j (i, j = 1, . . . , n) and i < 2n + 1 — j 
(i = 1, . . . , n, J = n + 1, . . . , 2n). We choose the value of each such , excluding those on the diagonal 
or anti-diagonal, according to the geometric distribution 

Pr(a;^j = fc) = (1 - qiqj){qiqjf 

On the anti-diagonal we choose 

V ( M / k even 

I 0, otherwise 

and on the diagonal 

Vv{x,,, = k) = [l - aq,){aq,f (i = l,...,n). (6.2) 

The bijection then implies that the probability X maps to a set of at most 2n non-intersecting u/rh 
lattice paths of final displacement 2A is equal to 

n n n+i 

n((l-9?)(l-«*) n (l-?.?.))"^" ^'4a (9i,---,'72„) (6.3) 

where 



2n 

s-,^-(gi,...,feO:=E*n'^l^'^'^''^'^ (6.4) 



with Xi{j) denoting the number of vertical steps at x = j — 1 contained in the level-^ path and the asterisk 
denoting that the sum is over all self dual u/rh non-intersecting lattice paths with final displacement 2 A. 
Because for such self dual lattice paths — Xi{2n + 1 — j) (recall the discussion below (|4.ip ) we have 



2n 



This allows us to set 



q2n+i~i (i = 1, • ■ • iJT-), 



and so with 



6.311 reads 



.l,^-(gi,...,g„) :=E*4f (gi,...,g2„) ^ l[lf'^' ''^'^ (6-5) 



9i=92Ti + l- 



- n (1 - 9.<?,)a^-^^-'^"'^^slf-(gi, . . . ,g„). (6.6) 

As noted in [2], the polynomial S2_^' is expressible in terms of Schur polynomials. To understand 
this point, one must first establish a relation between self dual lattice paths, represented as self dual 
tableaux, and domino tableaux. Regarding the latter, consider the diagram of a partition 2A. Define a 
domino tableau, of shape 2A with content from {n + 1, . . . , 2n}, as a tiling of the diagram by dominos with 
the dominos numbered from the set {71 + 1, ... , 2n} (each number repeated twice to fill the two squares 
of the domino) such that the number given to different dominos strictly increase down columns and 
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Figure 5: The self dual tableaux of shape 42 and content 4, the corresponding domino tableaux 
and the corresponding pair of tableaux of shape 2, content 2, and shape 1, content 1. 



weakly increase along rows. It is a known result (see e.g. [J^) that there is a bijection between self dual 
tableaux of shape 2A, content 2n, and domino tableaux of shape 2A with content from {n + 1, . . . , 2n}. 
In particular, to construct a domino tableau from a self dual tableau P say, one applies in succession 
the Schiitzenberger evacuation operation (see e.g. |20|). and the operation of removing the last square 
displaced in this operation. The domino formed by the evacuated and removed squares is numbered by 
the number of the removed square, which will be between n + 1 and 2n (note that P being self dual 
implies the sum of the entries of the evacuated and removed squares is 2n + 1). The procedure is repeated 
until all dominos have been identified and numbered (see Figure [S] for an example of the end product of 
this mapping). 

Having established the bijection between self dual tableaux and domino tableaux, one now makes 
use of a bijection between domino tableaux of shape 2A and content from {n + 1, . . . , 2n}, and pairs of 
semi-standard tableaux (P, Q) of shape (/i, k) each of content n with 

shapeP = (Ai, A3, . . . , Am) A+, shape Q = (A2, A4, A^^i) =: A^ (6.7) 

where m equals the length of 2A if the latter is odd, and one minus the length if it is even (and then 
Am = 0). To construct P (Q), remove all columns from the domino tableau for which the absolute value 
of the column number minus the length of the column is odd (even), and remove all even (odd) numbered 
rows. Finally subtract n from each of the numbers (an example of the result of this mapping is given 
in Figure [5]). Recalling the definition (|6.5p of S2x' and the definition (|2.4p of the Schur polynomials, it 
follows that 

S2\{Qi, ■■■,qn) = sx+{qi, ■ ■ .,qn)sx-{qi, . . .,?„). (6.8) 
Consequently the probability (j6.6p can be written as 

n n 

[|(1 - aq^) n (1 - 9^'Zj)a^^'=^'~'^'"'^^SA+ (gi, • ■ ■ , qn)sx- {qi, . . . , <7„). (6.9) 

In the recent work [9J (|6.9p was stated without derivation as being equal to the probability that the 
2n X 2n matrix X, symmetric about both the diagonal and anti-diagonal, and with elements distributed 
according to (|6.ip and (|6.2p , maps under the RSK correspondence to a set of at most 2n non- intersecting 
u/rh lattice paths of final displacement 2A. This is precisely the result derived here. 
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From the definition (|2.4p of the Schur polynomials in terms of non-intersecting lattice paths, it is easy 
to see that the well known identity 

X! a^J=i'"^'' ^'sx-{qi,---,q,i) ^ sx+{qi,---,qn,a) 

A+ fixed 

holds. Thus the marginal probability of A+ in (|6.9p is equal to 

n n 

Y[{i-aqi) Y[{^-qiq])sx+{qi,---,qn)sx+{qi,---,qn,a). (6.10) 

This in turn implies that with 

Lfn ■= max ^ 

(l,l)u/rh(2n,2n) 

we have 

n n 

Pr(Lf„ < 21) = Pr(if„ < 2/ + 1) = - aq.,) [] (1 - q,q,) 

1=1 ij = l 

X X! sx+iqi,---,qn)sx+{qi,---,qn,a). (6.11) 

A+:Ai</ 

The sum in (|6.1ip is a special case of that in (|2.1ip — thus replace n i-^ n + 1 in the latter and set 
Oi — bi — qi (i = 1, . . . , n), a„+i = 0, fe„+i = a. It therefore follows from p.3p that 

Pr(Lf„ < 20 = Pr(Lf„ < 2/ + 1) 

n n ^ n I 

= 11(1 - ag,) n (1 - q^q:){ 11(1 + ae^^--) I] IT |1 + 1^^'''\')um (^'^^^ 

1=1 ij = l fc=l j = l fc=l ^ ■* 

This result is the special case /3 = (the effect of setting /3 = is to constrain the elements on the 
anti-diagonal of X to be even) of a result first derived in [5] using methods of symmetric function theory 
to sum over a /3-generalization of (|6.3p . 

7 Matrices with a point reflection symmetry 

The point {n + 1/2, n + 1/2) is at the centre of the region 1 < x,y < 2n. A point {x, y) in this region 
reflected about this central point maps to the point (2n + 1 — x,2n + 1 — y). We thus say that a 2n x 2n 
matrix X has a point reflection symmetry (about the point {n + 1/2, n + 1/2)) if Xij — X2n+i-i,2n+i~j 
= 1, . . . ,n) or consequently ii X = (X^)"^. For a matrix with this symmetry we can take as the 
independent elements the triangular region below the anti-diagonal i < 2n + 1 — j = 1, . . . ,2n) 
together with the portion of the anti-diagonal i — 2n + I — j {i, j = 1, . . . ,n). 

We seek the constraint on the pairs of paths {Pi,P2), with both Pi and P2 of the same final 
displacements, which according to the RSK mapping are in correspondence with matrices X with 
a point reflection symmetry. We have already noted that with X mapping under RSK to (Pi,P2), 
X^ := [x2n+i-j,2n+i-i]i,]=i,...,2n maps to {P^^P") whilc X'^ := [xj,i]i^j=i^,,,^2n maps to {Pi , P^^) and 
hence matrices with the point reflection symmetry X ~ (X^)'^ map to a pair of u/rh lattice paths of the 
same final displacements constrained so that 

Pi = Pi^, P2=P^- (7.1) 
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We choose the independent entries of X according to the geometric distribution 

Pr(xjj = fc) = (1 - mj){q.iqj)^ 

where to be compatible with the point reflection symmetry we require g2n+i-i — Qi {i — 1, ■ ■ ■ ,n). With 
this specification it foUows from the bijection that the probabihty X maps to a pair of u/rh lattice paths 
of final displacement /i is equal to 

, 2 



Yl il-mj)sl^'^-{qi,...,qn) 



We know from (|6.8p that when the length of the parts of /i are all even, s^'^- can be given in terms 
of Schur polynomials. This is also true in the general case [H]. One again proceeds by noting that 
there is a bijection between a general self dual tableaux of shape A and domino tableaux. An immediate 
consequence is that unless A admits a domino tiling — for which the necessary and sufficient condition 
is that the number of points in the diagram of A with i + j even is equal to the number of points 
with i + j odd — one has s^'^' — 0. It is also true that domino tableaux are in bijective correspondence 
with pairs of semi-standard tableaux of shape (^'"^ /i'-^-'), l/x^"^! + l/i*-^-*! = each of content n where 
fi^^'' and fi^^^ are the so called 2-quotient of the partition fi. Regarding the latter, let fi he a, partition 
of length m. Add to fj, the partition Sm ■= ijn — 1, m — 2, . . . , 1, 0), and from this construct two new 
partitions fi^^"^ of lengths m^^\m,^^\ the first consisting of the even parts of fi + 6m, and the second 
the odd parts of /i + 6m- The 2-quotient is the pair of partitions /i^"-*, /i*-^^ specified by [17] 

^(0) = ^(o)/2 _ s^^,,, , = {f,W + l)/2 - ■ 

Analogous to (|6.8p one therefore has that if /i admits a domino tiling, then 

sll'^'iqi, . . . = s^(o)(gi, . . . ,g„)s^(i)(gi, . . . ,g„). (7.2) 

Regarding the converse of this statement, it's easy to see that the 2-quotient of a partition which admits 
a domino tiling is unique, while the 2-quotient of a partition which docs not admit a domino tiling 
coincides with the 2-quotient of a partition which does. Hence, given arbitrary partitions fi'--^\ /i*-^-* there 
is a unique partition /i which admits a domino tiling and is such that ()7.2|) is satisfied. Furthermore, from 



the definition of a 2-quotient, if /ui < 21 then < I and fj,'"^^ < I, while if /ii < 2/ + 1 then fi^^ < I + 1 
and /i^^'' < or ^J^-* < I + 1 and ^f''^ < I. Thus with 

Lfn max Xi^j 

(l,l)u/rh(2n,2n) 

we have 



n 

Pr(Li|, < 2/) = n (1 - l^^jf E • • ■ ' 9« 



n(l-*9,) E =(Pr(i°<0) ,,,,,,, (7-3) 

where the final equality follows upon comparison with (|2.1ip , and 



2 

Pr(Li}, < 2Z + 1) = n (1 - q,q,f ^ [sf{q,, . . . , g„)) 

n 

= W{^-mjf E is^{qi,...,qn)f Y {s^{qi,...,qn)) 

i-j — l k:k.i<1-\-1 k.:ki<1 



Pr(i°<; + 1) Pr(L,°<0 ^ ^ . (7.4) 
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The results (|7.3p and (|7.4p were stated without derivation in [3J . The derivation given here uses the strat- 
egy outUned in [5] to derive (|7.3p and (|7.4p in the exponential limit which corresponds to the appropriately 
symmetrized form of the Hammersley process. 
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